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Abstract 

An  algebraic  model  of  the  membrane  electrode  assembly  of  the  direct  methanol  fuel  cell  is  developed,  which  considers  the  simultaneous  liquid 
water  and  methanol  crossover  effects,  and  the  associated  electrochemical  reactions.  The  respective  anodic  and  cathodic  polarization  curves  can  be 
predicted  using  this  model.  Methanol  concentration  profile  and  flux  are  correlated  explicitly  with  the  operating  conditions  and  water  transport  rate. 
The  cathode  mixed  potential  effect  induced  by  the  methanol  crossover  is  included  and  the  subsequent  cell  voltage  loss  is  identified.  Water  crossover 
is  influenced  by  the  capillary  pressure  equilibrium  and  hydrophobic  property  within  the  cathode  gas  diffusion  layer.  The  model  can  be  used  to 
evaluate  the  cell  performance  at  various  working  parameters  such  as  membrane  thickness,  methanol  feed  concentration,  and  hydrophobicity  of  the 
cathode  gas  diffuser. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

One  critical  issue  of  the  research  on  direct  methanol  fuel 
cell  (DMFC)  is  to  reduce  methanol  crossover  through  the  poly¬ 
mer  electrolyte  membrane  [1 — 4],  The  permeation  of  methanol 
not  only  reduces  the  fuel  efficiency  but  also  causes  a  voltage 
decline  in  the  cathode  due  to  the  internal  short-circuit  formed 
by  the  parasitic  reactions  of  methanol  oxidation  and  the  oxygen 
reduction  [5-13].  Another  requirement,  especially  for  aportable 
DMFC,  is  the  retardation  of  water  transport  through  the  poly¬ 
mer  membrane.  Due  to  the  electro-osmotic  drag  of  water  across 
the  membrane,  the  use  of  concentrated  methanol  solution  in  the 
anode  compartment  is  not  feasible  [14-18].  Since  water  is  one  of 
the  reactants  for  methanol  oxidation,  the  anodic  limiting  species 
may  be  switched  from  methanol  to  water  once  water  loss  is 
severe  in  the  anode  compartment.  Condenser  is  needed  to  col¬ 
lect  water  in  the  cathode  and  to  recycle  it  to  the  anode;  thus,  the 
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design  of  the  portable  DMFC  becomes  much  complicated  and 
costly.  The  cathode  gas  diffusion  layer  is  normally  hydrophobic 
treated  to  induce  a  capillary  pressure  to  lessen  water  permeation 
from  the  anode.  An  extension  of  this  concept  is  to  fabricate  an 
extra  thin  micro-porous  layer  between  the  catalyst  layer  and 
the  cathode  gas  diffusion  layer  [16-20].  With  the  incorpora¬ 
tion  of  a  high  hydrophobic  and  low  porosity  carbon  layer,  the 
induced  capillary  pressure  is  further  extended,  while  the  impetus 
to  the  oxygen  diffusion  is  not  serious  due  to  the  thinness  of  the 
micro-layer. 

The  modeling  works  on  the  DMFC  with  attention  on  the 
methanol  crossover  can  be  found  in  [5-13] .  Scott  and  co-workers 
[6,7]  treated  the  cathode  voltage  loss  due  to  methanol  crossover 
as  a  fixed  value.  Kulikovsky  [8,9]  considered  the  effect  of 
channel  flow  on  the  cathode  mixed  potential,  and  found  that 
a  narrow  zone  of  local  current  bridge  which  short-circuits  the 
cell  at  a  small  current  load.  The  inclusion  of  water  transport 
phenomenon,  in  addition  to  the  methanol  crossover,  was  devel¬ 
oped  only  recently  [21-28].  It  is  noted  that  most  water  transport 
included  models  are  based  on  the  computational  fluid  dynamics 
except  a  few  analytical  works  in  Refs.  [21,22],  Sandhu  et  al. 
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Anode  Diffuser  Membrane  Cathode  Diffuser 


Anode  reaction :  CH3OH  +  H20  — »-  C02  +  61^+  6e' 

Cathode  reaction  :  02  +  4H*  +  4e' — ►  2HzO 

CH3OH  +  H20  -*■  C02  +  6H++  6e' 

Fig.  1 .  Schematic  of  the  membrane  electrode  assembly  of  DMFC. 


[21]  algebraically  formulated  the  methanol  and  water  crossover 
equations  but  only  limited  to  the  membrane  region  without  con¬ 
sideration  of  the  respective  anode  and  cathode  electrochemical 
reactions.  The  analytical  derivation  by  Shi  et  al.  [22]  was  also 
on  the  membrane  region  that  focused  on  the  water  flow  and 
membrane  dehydration;  methanol  crossover  and  electrochemi¬ 
cal  kinetics  were  not  included. 

The  purpose  of  present  work  is  to  develop  an  efficient  alge¬ 
braic  model  to  elucidate  the  simultaneous  water  and  methanol 
crossover  phenomena  as  well  as  their  influence  on  the  cell  per¬ 
formance.  To  limit  the  scope  of  this  study,  only  membrane 
electrode  assembly  (MEA)  is  considered,  while  the  flow  chan¬ 
nel  effect  is  not  included.  The  interplay  of  mass  transport  and 
associated  electrochemical  reactions  is  stressed.  Algebraic  while 
mechanistic -based  model  gives  an  explicit  relationship  between 
the  important  operating  parameters  and  the  physical  character¬ 
istics  of  the  MEA. 


2.  Mathematical  model 


2.1.  Liquid  water  flow 


Water  conservation  within  the  anode  diffuser  (AD): 


where  water  flux  is  governed  by  the  Darcy’s  law: 
k^dPi 


dx 


(1) 


(2) 


JV/AU  _ 

- - 


kp  D  is  the  liquid  permeability  in  the  anode  diffuser,  v  (=/r/ p)  the 
liquid  kinematic  viscosity  defined  as  the  ratio  of  liquid  viscosity 
and  molar  density,  Pi  is  the  local  liquid  pressure. 

Water  conservation  within  the  membrane  (M): 


-  =  0 


(3) 


dx 

The  water  flux  in  membrane  phase  is  driven  by  the  pressure 
gradient  and  the  electro-osmotic  drag: 


^ 


n?2U7 


(4) 


Water  conservation  within  the  cathode  diffuser  (CD): 


Water  flux  is  caused  by  hydraulic  pressure  gradient: 
»rCD  CCdP, 


(5) 


(6) 


*pCD  is  the  permeability  of  water  in  the  carbon  fiber  matrix. 
Within  the  porous  structure  of  the  cathode  diffusion  layer,  oxy¬ 
gen  gas  diffuses  in  the  non-flooding  fraction  of  the  void  space. 
The  permeability  of  water  transfer  is  modified  by  the  relative 
permeability,  k^D,  expressed  as  a  function  of  water  saturation 
=  .v3 4 5  [29] .  Water  saturation  s  is  determined  by  the  capillary 
equilibrium  between  the  gas  and  liquid  phases,  given  as  follows 
[29]: 


Schematic  of  the  studied  region  is  depicted  in  Fig.  1.  The 

assumptions  made  in  this  one-dimensional  algebraic  model  are 

listed: 

(1)  The  anode  and  cathode  catalyst  layers  are  treated  as  planes 
of  no  thickness.  That  is,  possible  concentration  variations 
within  the  catalyst  layers  are  not  considered. 

(2)  The  membrane  is  fully  hydrated. 

(3)  Electrochemical  kinetics  of  methanol  oxidation  and  oxygen 
reduction  is  described  by  the  Tafel  equation. 

(4)  Because  the  over-potential  for  the  methanol  oxidation  at 
cathode  is  high,  it  is  assumed  that  the  permeated  methanol 
is  consumed  instantaneously  at  the  membrane/cathode  dif¬ 
fuser  interface  [11,13]. 

(5)  The  effect  of  CO2  gas  bubbles  generated  in  the  anode  cata¬ 
lyst  layer  is  neglected. 


Pc  =  Pc  —  P\  =  0  cos  9C 


(7) 


Pc  is  the  capillary  pressure,  Pc  the  cathode  gas  pressure,  £cd  the 
porosity  within  cathode  porous  medium,  9C  the  contact  angel, 
o  the  surface  tension,  and  J(s)  is  the  Leverette  function  of  the 
water  saturation  [30]: 

J(s )  =  1 .4 17. v  —  2.12s1 2  +  1.263s3  for  hydrophobic  medium, 
9C  >  90°  (8a) 


J{s)  =  1.417(1  -  s)  -  2.12(1  -  s)2  +  1.263(1  -  s)3 

for  hydrophilic  medium  9C  <  90°  (8b) 

Water  balances  on  the  anode  reaction  plane  (AD/M  interface) 
and  the  cathode  reaction  plane  (M/CD  interface)  are  given, 
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respectively: 


juCD  _  nM  ,  7o 
ivH20  —  ^vHtO  +  777 


h2°  ^  2F 


/leak 

'  6  ~F 


(9) 


(10) 


lo  —  I  +  /leak 


/  70  X  _  pcosOc^^T 

[Nh2o  +  3f)~  (0.2415%/( 

— 0.6676.v'^cd  +  0.6315sm/cd) 

for0c  <  90° 


(15b) 


I  is  the  cell  current  density,  70  the  partial  current  density  of  oxy¬ 
gen  reduction,  and  7ieak  is  the  leaking  current  density  caused  by 
the  oxidation  of  permeated  methanol  on  the  cathode  catalyst. 

An  overall  current  conservation  can  be  expressed  in  Eq.  (11), 
which  states  that  the  oxygen  consumption  is  balanced  by  the 
oxidation  of  methanol  in  the  anode  plus  that  occurred  in  the 
cathode: 


2.2.  Methanol  transport 


Mass  conservation  of  methanol  within  the  anode  diffuser 
(AD): 


(11) 


dx 

Methanol  flux  is  contributed  by  diffusion  and  convection: 


(16) 


ndCm 


(17) 


Integration  of  Eq.  (2),  one  obtains: 

v<^AD 


pAD/M  _  pA  _  _ 


(12) 


Integrating  of  Eq.  (4)  with  the  aid  of  Eqs.  (9)  and  (12): 
oM/cd  _  pA  ,  vLmI  (1  h2cA 


Lm  Lad 


In  combination  of  Eqs.  (7)  and  (13),  one  has: 


vLmI  ( 1 

1  l  6  ' 


Cm  ,  Lad 


£  \0'5 

7(.vm/cd) 


D^u  is  the  void  fraction  corrected  diffusivity  using  Brugger- 
man’s  correlation  [31].  For  the  dilute  solution,  mole  fraction 
of  methanol  can  be  approximated  by  £"Cm,  where  is  the 
molecular  weight  of  water. 

Eq.  (17)  can  be  integrated  with  the  boundary  condition  Cm  = 
C £  atx  =  0: 


AC 


(13) 


r(<2D0  +  A/ad) 

,  w,  [r(<Do  +  /VAD)x 
+  Cm6XP[ - £)AD - 

Cm  at  AD/M  interface  is  simply: 

^AD/M  _  _ 


(18) 


Wigo  +  Ng>) 


|l  -  exp  |— 


Ko  +  O^ad 

dad 


(14) 


+C^  exp 


S"(N£>0  +  N™)Lad 
D 


(19) 


where  sM/cd  is  the  water  saturation  at  the  M/CD  interface.  Mass  conservation  of  methanol  within  the  membrane  phase  (M) 


With  the  aid  of  Eqs.  (7)— ( 1 1 ),  Eq.  (6)  can  be  integrated  over  the 
cathode  diffusion  layer  by  letting  water  saturation  zero  at  the 
interface  of  the  gas  diffuser  and  flow  channel: 


pcos6>c(£cD^pD)°' 


3F  )  v  Lqd 

+0.848s^/cd  —  0.6315sm/CD) 

for  0C  >90° 


(-0.35425^/cd 


dx 

Methanol  flux  is  expressed  as  follows: 
N™  =  -Dmd^m  +  ?"Cm(^H20  +  N™) 


(20) 


(21) 


Methanol  balance  on  the  anode  catalyst  layer  (AD/M  interface) 


(15a) 


and 


(22) 
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Integration  of  Eq.  (21)  with  the  aid  of  Eqs.  (19),  (20)  and  (22), 
one  obtains: 


-  1/6  F 


'  r(<20  +  ^mD-  IPF)  \  Wh?0  +  Aft D) 


1  —  exp 


+  exp 


+  N*D)LAD 

?"(<2D0  +  N^)Lad\ 
dad 


-  //6E 


exp 


r(<Do  +  AftD  -  1/3  F)  J 
'mffo  +  ^D-  I/3F)(x  -  Lad) 


(23) 


Due  to  the  high  cathode  over-potential  for  methanol  oxidation, 
methanol  is  consumed  instantaneously  at  the  M/CD  interface 
[11,13],  that  is,  Cm(jt  =  LAD+EM)  =  0,  we  have: 


0  = 


Aftu  -  I/6F 


U"(< 


?"(<o  + 

dad 


+  AftD-  I/3F) 

1  —  exp 

Aft13  -  7/6E 

”r(<2D0  +  /VnAD-//3C) 

+  -  I/3F)Lm\ 


Aftu 

o  +  Aft°) 


(24) 


2.3.  Oxygen  mass  transport 


Oxygen  conservation  within  the  cathode  diffuser  (CD): 


d  A£D 
dx 


=  0 


(25) 


Oxygen  flux  is  expressed  as  the  effective  diffusion  in  porous 
medium,  expressed  in  the  Bruggerman’s  relation: 


/V0CD  =  -A,[£cd(1  -  .s-m/cd)]1'5  (26) 

Although  the  void  space  available  for  the  gas  passage  is  a  func¬ 
tion  of  local  water  saturation  s,  we  adopt  the  smallest  void 
fraction  £cd(  1  —  ^m/cd)  as  a  first  approximation.  Thus,  the  resis¬ 
tance  to  oxygen  diffusion  by  water  flooding  is  over-estimated  in 
the  present  study. 

<D  is  related  to  the  oxygen  partial  current  density: 


/V0CD 


Io 

4F 


(27) 


Oxygen  concentration  at  M/CD  interface  is  derived  in  combina¬ 
tion  of  Eqs.  (25)-(27): 

rM/CD  _  rb _ IqLcd _ 

0  0  4FD0[£Cd(1— sm/cd)]15 

2.4.  Electrode  kinetics  expressions 


The  oxygen  partial  current  density  at  cathode  is  expressed  in 
Tafel  equation: 


KCoM 


(29) 


K  is  the  solubility  constant  of  oxygen  between  gas  and  aque¬ 
ous  phases  (1  —  sm/cd)  accounts  for  the  fraction  of  available  gas 
contact  area.  /00,ref  =  *oo,ref  fwc«cs  in  which  i00,ref  is  the  reference 
exchange  current  density  of  oxygen  reduction,  me  the  cathode 
catalyst  loading,  and  ac  is  the  active  surface  area  per  unit  cata¬ 
lyst  mass.  The  anode  methanol  oxidation  current  density  can  be 
expressed  similarly: 


^AD/M 

I  =  7om,ref  ^  CXp 

Cm, ref 


(t£) 


(30) 


7om,ref  =  *om,ref  ^AflA,  in  which  i0m,ref  is  the  reference  exchange 
current  density  of  methanol  oxidation,  m\  the  anode  catalyst 
loading,  and  a\  is  the  active  surface  area  per  unit  mass  of  the 
anode  catalyst. 


2.5.  Cell  polarization  expression 


The  cell  voltage  is  calculated  by  subtracting  the  voltage  losses 
of  cathode  \rjc\,  anode  rja,  and  membrane  Lyjhcm  from  the  ther¬ 
modynamic  cell  potential  UQ  —  Um: 

Vcell  =  UQ-Um-  he  I  -  -  —  (31) 

Km 


2.6.  Solution  procedure 

The  cell  polarization  curve  is  calculated  in  a  potentiostatic 
mode.  For  a  given  Vceu,  variables  rjc,  0  a,  4  A^q,  sm/cd,  and 
Aft°  are  to  be  determined  by  Eqs.  (29)-(31),  (14),  (15),  and 
(24).  Note  that  C^CD  and  in  Eqs.  (29)-(30)  are  given 

exphcitly  by  Eqs.  (19)  and  (28),  while  /<,  in  Eqs.  (15),  (28)  and 
(29)  is  correlated  by  mass  conservation  over  the  entire  MEA: 

I0  =  6  Fn£°  (32) 

Newton-Raphson  method  [32]  is  used  to  solve  the  set  of 
equations  iteratively.  The  fixed  parameters  used  in  the  model 
are  fisted  in  Table  1  unless  otherwise  specified. 

3.  Results  and  discussion 

3.1.  Membrane  thickness  effect 

Fig.  2  depicts  the  cell  polarization  behavior  of  varied  mem¬ 
brane  thickness.  Accompanied  in  this  figure  is  the  corresponding 
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Table  1 

Fixed  mass  transport  and  kinetic  parameters 


Parameter 

Expression 

References 

OA 

52  x  104  cm2  g-1 

[33] 

ac 

52  x  104  cm2  g-1 

[11] 

CH+ 

KmRT(1000) 

[31] 

Cm, ref 

0.5  x  10-3  mol  cm-3 

[11] 

C0,ref 

0.496  x  10-6  mol  cm- 3 

Corresponding  to  the  solubility  of  1  atm  O2  at  353  K 

D0 

1  823(PcaTcb)°'333(Tcarcb)0-4167  yfiz+b 

[34] 

% 

1  x  10-4  cm2  s-1 

[35] 

d£d 

0.1274  x  10“4  cm2  s"1 

[11],  void  fraction  corrected  diffusivity 

1.2  x  10“6  A  cm-2 

[38] 

inm/p.f 

0.7630  x  10-7  Acm-2 

[11] 

«oo  ref 

0.3189  x  10-7  Acm-2 

[36] 

if 

6.875  x  10“9  cm2 

[29] 

6.875  xl0“9  cm2 

Chosen  the  same  as 

k™ 

2  x  10-15  cm2 

An  average  of  Refs.  [28,37] 

K 

W  exp  (=^  +  14.1) 

[33] 

Lad 

0.026  cm 

[31] 

Lq  d 

0.026  cm 

[31] 

Cm 

0.0051,  0.0127,  0.023  cm  for  Nafion  112,  115,  and  117 

[31,36] 

„H  2° 

2.5 

[38] 

Ma 

32  g  mol-1 

[34] 

Mb 

28  g  mol-1 

[34] 

Pc  a 

49.7  atm 

[34] 

Pcb 

33.5  atm 

[34] 

Tea 

154.4  atm 

[34] 

Ub 

126.2  atm 

[34] 

t/mo 

0.03  V 

[13] 

um 

Cmo  +  F  In  (  (Cn^reflOOO)1/6  ) 

[31],  (reference  pressure  of  CO2  =  1  atm) 

Uoo 

1.23  -  0.9  X  10“3  (T  —  298) 

[31] 

Uo 

U00  +  In  CH+ 

[31],  (reference  pressure  of  O2  =  1  atm) 

t*a 

0.5 

[27,28] 

1 

[28,33] 

Km 

0.6842  x  10_1  S  cm-1 

[11,31,39] 

ecD 

0.6 

[23] 

a 

62.5  gs-a 

[29] 

M 

3.565  x  10“3  gem”1  s”1 

[29,31] 

P 

0.054  mol  cm-3 

[31] 

K" 

18gmol_1 

Water  molecular  weight 

methanol  crossover  rate  expressed  as  the  leaking  current  density. 
It  is  seen  that  methanol  permeation  is  significant  for  the  thinner 
membrane,  which  is  declined  gradually  when  the  cell  limiting 
current  density  is  approached.  Usually  highest  leaking  current 
density  is  observed  at  the  open  circuit  voltage;  however,  there 
are  occasions  that  a  maximum  methanol  crossover  rate  appears 
in  the  polarized  region.  We  will  come  back  to  this  point  later 
when  we  discuss  the  methanol  feed  concentration  effect.  The 
coincidence  of  the  vanishing  in  leaking  current  density  with 
the  arrival  of  cell  limiting  current  density  implies  the  domi¬ 
nance  of  anodic  process  for  the  prescribed  parameters.  There 
is  a  decline  of  cell  voltage  near  OCV,  covering  the  activation 
till  the  emerging  of  ohmic  resistance  controlled  regime,  when 
a  thin  membrane  is  used.  The  voltage  loss  of  thin  membrane  is 
originated  from  the  mixed  potential  effect.  The  higher  methanol 


crossover  rate  induces  an  internal  short-circuit  on  the  cathode 
electrode;  thus,  the  available  cathode  potential  is  partially  con¬ 
sumed  by  the  oxidation  of  crossed  methanol  on  the  cathode.  As 
the  discharge  rate  increases,  the  crossover  effect  is  diminished 
while  the  membrane  ohmic  resistance  dominates  as  more  sig¬ 
nificant  for  the  thicker  membrane.  For  a  better  illustration  of  the 
relative  contributions  of  anode,  cathode,  and  ohmic  resistance, 
Fig.  3  depicts  the  cumulative  over-potentials  for  Nafion  112  and 
1 17.  As  shown  in  this  plot,  Nafion  1 12  has  larger  cathode  over¬ 
potential  \rjc\  than  that  of  1 17  due  to  the  mixed  potential  effect; 
the  difference  decreases  at  larger  discharge  rate.  It  is  interesting 
to  note  that  when  the  anodic  over-potential  rjn  is  included,  that  is, 
the  portion  accounted  by  U0  —  Um  —  \0c\  —  Va>  Nafion  112  has 
an  overall  higher  over-potential  than  Nafion  1 17.  This  is  caused 
by  the  always  less  methanol  concentration  available  on  the  anode 
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lcen  /A  cm"2 

Fig.  2.  Cell  polarization  curve  with  the  associated  leaking  current  density  for 
different  membrane  thickness.  Cathode  feed:  air,  anode  feed:  methanol  0.5  M, 
cathode  gas  diffuser  contact  angle:  120°,  PA  =  PC  =  1  atm,  T  =  353  K. 

catalyst  when  the  thinner  membrane  employed.  Lower  concen¬ 
tration  induces  higher  anode  electrode  over-potential  rjd.  This 
point  will  become  clear  when  methanol  concentration  profile  is 
shown  later  in  Fig.  4.  Finally,  as  the  ohmic  potential  Tm/cen//cm 
is  included,  the  thicker  membrane  has  the  largest  polarization  in 
the  region  of  high  discharge  rates. 


Fig.  3.  Cumulative  over- voltages  for  Nations  112  and  117.  Cathode  feed: 
air,  anode  feed:  methanol  0.5 M,  cathode  gas  diffuser  contact  angle:  120°, 
pA=pC  =  latmj  7=353 k. 


Nation  117 


Anode  diffuser  Nation  1 1 2 

< - X - > 


Fig.  4.  Methanol  concentration  profiles  within  the  anode  diffuser  and  the  mem¬ 
brane  for  Nations  112  and  1 17.  Cell  voltages  are  specified  at  0.6,  0.4,  and  0.2  V. 
Cathode  feed:  air,  anode  feed:  methanol  0.5  M,  cathode  gas  diffuser  contact 
angle:  120°,  PA  =PC  =  1  atm,  F=  353  K. 

Detailed  methanol  concentration  profiles  are  given  at  three 
cell  voltages  depicted  in  Fig.  4.  It  is  shown  at  a  specified 
potential,  thicker  membrane  preserves  a  higher  methanol  con¬ 
centration  at  AD/M  interface.  The  less  available  concentration 
for  thinner  membrane  confirms  the  observation  of  higher  anode 
over-potential  r]d  in  Fig.  3.  It  is  the  larger  mass  diffusion  of 
methanol  within  thinner  membrane  decreases  methanol  concen¬ 
tration  on  the  anode  catalyst  and,  in  turn  increases  rj a.  Zhang 
and  Wang  [10]  defined  “the  crossover  over-potential”  as  “the 
difference  between  the  cathode  over-potential  of  a  real  DMFC 
and  the  ideal  case  with  no  crossover”.  Although  the  crossover 
phenomenon  can  be  characterized  by  the  rise  of  cathode  over¬ 
potential,  it  is  shown  in  the  present  study  that  crossover  affects 
the  anode  potential  as  well,  and  should  be  noticed. 

Water  saturation  profiles  across  the  cathode  gas  diffuser 
in  Fig.  5  are  specified  at  the  cell  voltages  of  0.6,  0.4  and 
0.2  V.  Zero  saturation  is  presumed  at  the  diffuser/flow  channel 
interface  in  the  calculation.  Capillary  pressure  equilibrium  and 
water  transport  determines  the  saturation  profiles.  As  expect,  a 
larger  over-voltage  induces  higher  water  saturation  caused  by 
the  significant  electro-osmotic  water  drag  through  membrane 
and  water  generation  at  the  cathode  catalyst  layer.  Certain  fea¬ 
ture  is  observed  in  the  water  saturation  profiles.  That  is,  near 
OCV  (0.6  V),  thicker  membrane  induces  higher  water  satura¬ 
tion  within  the  cathode  diffuser.  This  trend  reverses  at  high 
over-voltage  (0.2  V).  Such  behavior  is  in  parallel  with  the  cell 
polarization  curves  in  Fig.  2.  For  a  specified  voltage,  the  current 
response  is  in  close  association  with  the  water  hold  up  in  the 
cathode  diffuser.  In  addition  to  the  water  generation  of  oxygen 
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Fig.  5.  Water  saturation  profiles  within  the  cathode  diffusion  layer  (CD)  with 
varied  membrane  thickness.  Cell  voltages  are  specified  at  0.6,  0.4,  and  0.2  V. 
Cathode  feed:  air,  anode  feed:  methanol  0.5  M,  cathode  gas  diffuser  contact 
angle:  120°,  PA  =P°  =  1  atm,  T=  353  K. 

reduction  in  the  cathode  catalyst  layer,  Eq.  (4)  indicates  that  the 
electro-osmotic  drag  of  water  in  membrane  is  in  proportion  to 
the  cell  discharge  rate. 

3.2.  Effect  of  methanol  concentration 

Fig.  6  displays  the  cell  performance  at  different  methanol 
feed  concentrations.  The  coincidence  of  the  approaching  of  cell 
limiting  current  density  with  the  vanishing  of  methanol  leaking 
current  density  suggests  the  dominance  of  the  anodic  process 
for  the  given  concentrations  and  operating  conditions.  As  expect, 
higher  cell  current  density  occurs  when  a  high  methanol  concen¬ 
tration  provided.  The  limiting  cell  current  density  is  proportional 
to  the  methanol  concentration  used.  However,  there  is  peculiar 
feature  appears  when  various  feed  concentrations  are  compared. 
There  is  an  enhanced  proportionality  with  methanol  concentra¬ 
tion  as  regard  to  the  limiting  current  density.  For  instance,  the 
limiting  current  density  is  over  0.6  A  cm-2  at  1  M  methanol  fed, 
which  is  more  than  twice  of  that  at  0.5  M.  The  limiting  current 
densities  when  convective  mass  transport  excluded  can  be  cal¬ 
culated  according  to  6/7D^d(C^1/ Lad),  which  give  the  values 
of  0.06514,  0.2606  and  0.5211  A  cm-2  for  methanol  concen¬ 
trations  0.125,  0.5  and  1 M,  respectively.  The  model  calculated 
limiting  current  densities  when  convective  mass  transport  con¬ 
sidered  in  Fig.  6  are  0.06639, 0.2824  and  0.6218  A  cm-2  for  the 
given  concentrations.  The  enhancement  is  a  direct  consequence 
of  the  water  flow-induced  convective  mass  transfer  on  methanol, 
as  described  in  Eqs.  (17)  and  (21).  Fig.  6  also  shows  that  higher 
methanol  concentration  induces  a  higher  leaking  current  density. 
Although  not  clear  in  this  figure,  there  are  interesting  maximum 
leaking  current  densities  appearing  at  0.46  and  0.24  V  for  the 
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Fig.  6.  Cell  polarization  curves  with  the  corresponding  methanol  leaking  current 
densities  at  various  methanol  feed  concentrations.  Cathode  feed:  air,  membrane 
Nation  117,  cathode  gas  diffuser  contact  angle:  120 °,PA=PC  =  1  atm,  7=  353  K. 

methanol  feeds  of  0.5  and  1  M,  respectively.  Fig.  7  gives  a  mag¬ 
nified  scale  about  the  leaking  current  density  as  a  function  of  cell 
current  density.  There  is  a  initial  rise  of  methanol  crossover  with 
the  cell  current  density  when  the  higher  methanol  concentration 
is  used.  The  nonlinear  behavior  has  to  do  with  the  increase  of 
convective  contribution  of  methanol  transport  as  suggested  in 


Fig.  7.  Leaking  current  density  as  a  function  of  cell  current  density  at  various 
methanol  concentrations.  Cathode  feed:  air,  membrane  Nafion  1 17,  cathode  gas 
diffuser  contact  angle:  120°,  PA=PC  =  1  atm,  T=  353  K. 
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Fig.  8.  Cell  voltage  vs.  water  flux  rates  within  anode  diffuser  (AD),  mem¬ 
brane  (M),  and  cathode  diffuser  (CD)  at  various  feed  concentrations.  Cathode 
feed:  air,  membrane  Nation  117,  cathode  gas  diffuser  contact  angle:  120°, 
pA=pc  =  latm,  T=353K. 

the  second  term  on  the  right  hand  side  of  Eq.  (21).  Although  the 
diffusive  mass  transfer  decreases  steadily  with  7cen,  there  could 
be  a  situation  where  the  enhanced  methanol  convection  by  water 
flow  surpasses  the  decline  of  diffusion  so  that  an  overall  increase 
of  beak  occurs.  Eq.  (4)  expresses  the  linear  increase  of  water  flux 
with  the  cell  current  density.  At  a  later  stage,  both  mechanisms 
of  methanol  transport  decrease  due  to  the  consumption  of  reac¬ 
tant  at  the  anode  catalyst.  The  occurring  of  methanol  crossover 
maximum  has  to  do  with  the  anode  catalyst  activity,  methanol 
feed  concentration,  and  the  transport  parameters  within  the  diffu¬ 
sion  media.  Similar  observation  was  found  by  other  researchers 
[27,28],  and  even  in  the  direct  ethanol  fuel  cell  systems  [40,41], 

Fig.  8  illustrates  the  respective  water  fluxes  through  the  anode 
diffuser  (AD),  membrane  (M),  and  cathode  diffuser  (CD)  with 
the  cell  voltages.  The  water  discharge  rate  is  in  parallel  with 
the  polarization  curve.  Eqs.  (2),  (9)  and  (10)  determine  the  dis¬ 
tribution  of  water  fluxes  within  various  media  of  the  MEA. 
Though  not  shown  here,  higher  methanol  concentration  also 
induces  higher  water  saturation  in  cathode  diffuser  and  higher 
back-pressure  through  the  membrane. 

3.3.  Inference  of  contact  angle  in  the  cathode  diffuser 

In  addition  to  the  permeability  and  porosity  of  the  gas  diffu¬ 
sion  layer,  the  water  transport  depends  on  the  hydrophobicity  of 
the  layer,  which  is  normally  characterized  by  the  contact  angle 
between  liquid  water  and  carbon  fiber  matrix  as  shown  in  Eq. 
(7).  Lim  and  Wang  [  1 8]  in  the  study  of  proton  exchange  fuel  cells 
showed  that  for  various  contents  of  fluorinated  ethylene  propy¬ 
lene  (FEP)  added  in  the  layer,  the  contact  angles  measured  were 
decreased  rapidly  with  temperature  showing  a  transition  from 


Fig.  9.  Cell  polarization  curves  with  the  corresponding  leaking  current  densities 
at  various  cathode  diffuser  (CD)  contact  angles.  Cathode  feed:  air,  membrane 
Nation  117,  methanol  0.5  M,  PA  =  Pc  =  1  atm,  T=  353  K. 

hydrophobic  ( 0C  >  90°)  to  hydrophilic  property  (9C  <  90°).  That 
is,  hydrophobic  property  of  the  diffuser  may  change  with  the 
operating  conditions.  Fig.  9  depicts  the  simulated  polarization 
curves  for  a  range  of  contact  angles  in  the  CD.  The  corresponding 
leaking  current  densities  are  included  as  well.  There  is  hardly  any 
significant  variation  of  cell  current  density  with  the  contact  angle 
under  the  specified  MEA  structure  and  mass  transport  parame¬ 
ters.  The  reason  for  that  is  the  predicted  saturation  within  the  CD 
is  low  (always  less  than  10%)  so  that  there  is  only  minor  pore 
blocking  effect  on  gaseous  oxygen  transport  in  the  CD.  In  addi¬ 
tion,  the  cell  discharge  is  mainly  limited  by  the  anodic  reaction 
so  that  characteristics  in  the  cathode  diffuser  have  insignificant 
influence  on  the  MEA  performance. 

Liquid  pressure  at  M/CD  interface  pM/CD  along  with  the  cell 
voltage  at  different  CD  contact  angles  is  demonstrated  in  Fig.  10. 
Liquid  pressure  at  AD/M  interface  is  also  included,  which  is 
nearly  constant  during  the  discharge  process.  That  is,  the  sec¬ 
ond  term  on  the  right  side  of  Eq.  (12)  is  negligible.  pM/CD  at 
OCV  depends  on  the  water  generation  at  mixed  potential  of  the 
cathode  and  the  associated  capillary  pressure  equilibrium  on  the 
water  saturation.  Positive  capillary  pressure  is  developed  for  the 
hydrophobic  medium,  while  negative  capillary  pressure  is  for 
the  hydrophilic  diffuser.  It  seems  the  magnitude  of  the  capillary 
pressure,  \PC\,  for  the  hydrophilic  medium  is  much  larger  than 
that  of  hydrophobic  medium.  In  addition,  the  liquid  pressure 
pM/CD  ajwayS  increases  with  the  cell  voltage,  as  the  result  of 
increased  water  saturation,  indifferent  to  the  hydrophobicity  of 
the  CD.  Thus,  a  positive  increase  of  (PM/CD  -  PAD/M)/Lm  with 
over-voltage  is  developed  for  the  hydrophobic  diffuser,  and  an 
increased  back-pressure  is  induced  as  indicated  in  the  first  term 
on  the  right  side  of  Eq.  (4).  The  corresponding  water  fluxes  in 
different  domains  of  the  MEA  are  demonstrated  in  Fig.  1 1 .  It  is 
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Fig.  10.  Liquid  pressure  at  membrane/cathode  diffuser  interface  vs.  cell  voltage 
at  various  cathode  diffuser  contact  angles.  Cathode  feed:  air,  membrane  Nation 
117,  methanol  0.5  M,  PA  =  P0  =  1  atm,  T=  353  K. 


clear  that  higher  water  flux  occurs  in  each  domain  when  the  CD 
becomes  more  hydrophilic.  As  discussed  earlier,  water  flood¬ 
ing  in  the  CD  is  not  severe  so  that  no  significant  impact  on  the 
polarization  curve  is  observed  in  Fig.  9.  However,  it  should  be 
mentioned  that  there  is  minor  difference  in  Fig.  9  at  high  dis¬ 
charge  rate  if  magnified  scale  is  plotted,  due  to  the  increased 


NH2q/  mol  cm'2  s'1 

Fig.  11.  Cell  voltage  vs.  water  flux  rates  within  anode  diffuser  (AD),  mem¬ 
brane  (M),  and  cathode  diffuser  (CD)  at  various  cathode  diffuser  contact  angles. 
Cathode  feed:  air,  membrane  Nation  117,  methanol  0.5  M,  PA  =  Pc  =  1  atm, 
T=353  K. 


back-pressure  within  the  membrane  as  justified  in  Fig.  10  and 
Eq.  (4).  Fig.  1 1  indicates  that  the  hydrophobic  layer  indeed  has 
a  better  capability  in  inhibiting  the  water  flow.  To  extend  the 
back-pressure  for  preventing  water  flux,  an  extra  thin  hydropho¬ 
bic  micro-layer  may  be  added  to  the  M/CD  interface  as  did  by 
Xu  et  al.  [20]  and  Pasaogullari  et  al.  [37]. 

4.  Conclusion 

Analytical  mass  transfer  equations  on  the  membrane  elec¬ 
trode  assembly  of  direct  methanol  fuel  cell  are  formulated,  which 
give  explicitly  the  inter-relationship  of  methanol  crossover 
and  water  transport  rate.  The  cell  voltage  drop  caused  by  the 
methanol  crossover  is  evaluated  by  the  mixed  potential  theory  on 
the  parasitic  oxygen  reduction  and  methanol  oxidation  occurring 
on  the  cathode  catalyst.  Water  flux  rate,  methanol  concentration, 
oxygen  gas  concentration,  liquid  pressure  distribution  and  water 
saturation  in  the  MEA  can  be  simulated  in  corporation  with  the 
given  electrochemical  kinetics. 

According  to  the  model  prediction,  larger  voltage  drop  occurs 
near  OCV  due  to  the  mixed  potential  effect  when  the  thinner 
membrane  is  used;  however,  the  voltage  gained  at  larger  polar¬ 
ization  benefits  from  the  smaller  ohmic  resistance  as  compared 
to  the  thicker  membrane.  It  is  also  pointed  out  the  methanol 
crossover  not  only  affects  the  cathode  over-potential  but  also 
varies  the  methanol  concentration  distribution  near  the  anode, 
and  alters  the  anode  over-potential.  In  the  study  of  methanol  feed 
concentration  effect,  it  is  shown  that  the  convective  contribution 
on  methanol  crossover  by  water  flux  may  surpass  the  decrease  of 
diffusion  mechanism  in  the  early  stage  of  cell  discharge,  caus¬ 
ing  a  maximum  on  the  leaking  current  density.  Hydrophobic 
treatment  on  the  cathode  diffusion  layer  seems  not  affect  the 
polarization  behavior  when  the  CD  water  saturation  is  not  high 
enough  to  block  the  gas  passage.  However,  the  water  flux  pre¬ 
dicted  is  indeed  a  function  of  surface  contact  angle,  and  less 
water  flux  is  observed  when  using  the  hydrophobic  layer  because 
of  the  induced  back-pressure  across  the  membrane. 
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Appendix  A.  Nomenclature 

a  a  active  area  per  unit  mass  of  anode  catalyst  (cm2  g-1) 

ac  active  area  per  unit  mass  of  cathode  catalyst  (cm2  g~ 1 ) 

CH+  proton  concentration  in  membrane  phase  (mol  dm-3) 

Cm  methanol  concentration  (mol  cm-3) 

Cm,ref  reference  methanol  concentration  (mol  cm-3) 

C0  oxygen  concentration  (mol  cm-3) 

C0, ref  reference  oxygen  concentration  (mol  cm-3) 
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Cj^  methanol  concentration  in  the  anode  flow  channel 
(mol  cm-3) 

C,r/M  methanol  concentration  at  anode  diffuser/membrane 
interface  (mol  cm-3) 

Cq  oxygen  gas  concentration  in  the  cathode  flow  channel 
(mol  cm-3) 

£,m/cd  OXygen  gas  concentration  at  membrane/cathode  dif¬ 


fuser  interface  (mol  cm-3) 

Dh+  proton  diffusivity  in  membrane  (cm2  s-1) 

Da  oxygen  gas  diffusivity  in  cathode  gas  diffuser  (cm2  s- 1 ) 
Dad  methanol  diffusivity  in  the  anode  diffuser  (cm2  s-1) 
methanol  diffusivity  in  the  membrane  (cm2  s-1) 

F  Faraday  constant  (96,500  C  mol- 1 ) 


?om,ref  reference  methanol  oxidation  exchange  current  density 
(A  cm-2) 

ioo,ref  reference  oxygen  reduction  exchange  current  density 
(A  cm-2) 

I  cell  current  density  (A  cm-2) 

/leak  crossover  current  density  (A  cm-2) 

70  oxygen  reduction  current  density  (A  cm-2) 

/om,ref  ?om,ref  (A  Cm-2) 

/oo,ref  Wef'AAc  (A  cm-2) 

kpD  anode  diffuser  liquid  permeability  (cm2) 
kp  0  cathode  diffuser  liquid  permeability  (cm2) 
membrane  liquid  permeability  (cm2) 
liquid  phase  relative  permeability  in  cathode 
diffuser 

K  solubility  constant  of  oxygen  in  liquid  phase 
Lad  anode  backing  layer  thickness  (cm) 

Lcd  cathode  gas  diffuser  thickness  (cm) 

Lm  membrane  thickness  (cm) 

m\  anode  catalyst  loading  (g  cm-2) 

me  cathode  catalyst  loading  (g  cm-2) 

Ma  molecular  weight  of  O2  (g  cm-2) 

Mb  molecular  weight  of  N2  (g  cm-2) 

rtj2°  electro-osmotic  drag  coefficient 

<D0  water  flux  in  anode  diffuser  (mol  cm-2  s-1) 

/V™0  water  flux  in  cathode  diffuser  (mol  cm-2  s-1) 

Ko  water  flux  in  membrane  (mol  cm  2  s  1 ) 

(VAD  methanol  flux  in  anode  diffuser  (mol  cm-2  s-1) 
methanol  flux  in  membrane  (mol  cm-2  s-1) 

A0cd  oxygen  gas  flux  in  cathode  diffuser  (mol  cm  2  s  1 ) 

Pc  capillary  pressure  (g  cm-1  s-2) 

Pca  critical  pressure  of  O2  (atm) 

Pcb  critical  pressure  of  N2  (atm) 

Pi  liquid  pressure  (g  cm-1  s-2) 

PA  anode  pressure  (g  cm-1  s-2) 

pAD/M  Hquid 

pressure  at  anode  diffuser/membrane 
(gcm-1s-2) 

Pc  cathode  pressure  (gem-1  s-2) 

pM/CD  liquid 

pressure  at  membrane/cathode  diffuser 
(gem-1  s-2) 

R  universal  gas  constant  (8.314  Jmol-1  K-1) 

Pg  universal  gas  constant  (82.06  atm  cm3  mol-1  K-1) 

,v  water  saturation  in  cathode  diffuser 


sm/cd  water  saturation  at  membrane/cathode  diffuser  inter¬ 
face 

T  absolute  temperature  (K) 

Pca  critical  temperature  of  O2  (K) 

Tc b  critical  temperature  of  N2  (K) 

f/m  reference  methanol  oxidation  open  circuit  potential  (V) 
Umo  standard  potential  of  methanol  oxidation  (V) 

U0  reference  oxygen  reduction  open  circuit  potential  (V) 
UQO  standard  potential  of  oxygen  reduction  (V) 

Vceii  cell  voltage  (V) 

Greek  letters 

aa  anodic  transfer  coefficient  of  methanol  oxidation 
ac  cathodic  transfer  coefficient  of  oxygen  reduction 
scd  void  fraction  in  cathode  gas  diffuser 

f"  molecular  weight  of  water  (g  mol- 1 ) 
r]a  electrode  over-potential  in  anode  (V) 

77c  electrode  over-potential  in  cathode  (V) 

9C  contact  angle  (rad) 

Km  proton  conductivity  in  membrane  phase  (S  cm-1) 
fx  viscosity  of  water  (g  cm-1  s-1) 
v  kinematic  viscosity  (g  cm2  s-1  mol-1) 

p  liquid  water  molar  density  (mol  cm-3) 

a  surface  tension  (g  s-2) 
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